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Abstract 

A systematic procedure to study one-dimensional Schrodinger equation 
with a position-dependent effective mass (PDEM) in the kinetic energy op- 
erator is explored. The conventional free-particle problem reveals a new and 
interesting situation in that, in the presence of a mass background, forma- 
tion of bound states is signalled. We also discuss coordinate-transformed, 
constant-mass Schrodinger equation, its matching with the PDEM form and 
the consequent decoupling of the ambiguity parameters. This provides a uni- 
fied approach to many exact results known in the literature, as well as to a 
lot of new ones. 
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1 Introduction 



In recent times, the concept of a position-dependent-effective-mass (PDEM) quan- 
tum Hamiltonian is rapidly gaining acceptance because of its increasing relevance 
in describing the motion of electrons in problems of compositionally graded crys- 
tals [1] (following our ability to fabricate semiconductor nanostructures), quantum 
dots [2], liquid crystals [3], etc. The appearance of PDEM is also well known 
in the energy density functional approach to the nuclear many-body problem [4] 
and its applications [5, 6] in the context of nonlocal terms of the accompany- 
ing potential. Other theoretical considerations where PDEM has been exploited 
include the derivation [7] of the underlying electron Hamiltonian from instanta- 
neous Galilean invariance and implementation of the path integral techniques [8] 
to calculate the Green's function [9] for step and rectangular-barrier potentials and 
masses. Further, PDEM has proved to be appealing in the construction of accept- 
able quantum mechanical systems by seeking exact solutions of the Schrodinger 
equation [10, 11, 12, 13, 14, 15, 16, 17, 18] by extending the already existing meth- 
ods of spectrum generating or potential algebras [19] and those of supersymmetric 
quantum mechanics [20, 21, 22, 23]. 

In this Letter we consider the PDEM formalism to address the issue of the 
generation of the corresponding potentials. We use the von Roos effective-mass 
kinetic energy operator [24], which has the advantage of an inbuilt Hermiticity, and 
we develop a PDEM scheme in very general terms that not only encompasses some 
of the previous ones but also brings out additional features not reported before. 

First, we consider the problem of a free particle in the framework of a first-order 
intertwining relationship and show that in a suitable mass background it acts as a 
bound system. 

Second, we observe that a Schrodinger equation with constant mass and having 
for its potential U can always be coordinate transformed in terms of some given 
mass function such that the resulting form matches with the one implied by PDEM. 
In this way the starting potential V of the PDEM scheme becomes expressible in 
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terms of the coordinate-transformed U and the one induced by the effective kinetic 
energy operator T. The latter is found to depend on two types of parameters: those 
entering the definition of the mass function and those coming from the ordering 
ambiguity of the momentum and mass operators. 



tive potential which a particle with a position-dependent mass experiences, it is not 
at all apparent how V stands out against the interplay of these ambiguity parame- 
ters vis-a-vis the ones of the coordinate-transformed U. In fact, as we shall see, in 
a PDEM problem the ambiguity parameters get decoupled in a natural way. This 
observation is new and of significance in that it allows V to become entirely identi- 
fiable with the coordinate-transformed U. The latter in the constant-mass limiting 
case goes over to a known form in the sense that it is either exactly solvable (ES) or 
quasi-exactly solvable (QES) or else conditionally exactly solvable (CES). We also 
discuss some consequences of our results. 



2 General Strategy and the Free-particle Prob- 



The most commonly used effective-mass kinetic energy operator is the two- 
parameter form given by von Roos [24], which in one dimension reads 



where p represents the momentum operator, m(x) is the position-dependent mass 
and the parameters a, (3, 7 are subject to the condition a + (3 + 7 = —1. Because of 
the noncommutativity of the momentum and (position-dependent) mass operators, 
maintaining Hermiticity of the kinetic energy operator is not trivial. In (1), T has 
been specifically designed to be Hermitian but is by no means unique. However, 
it turns out that other plausible forms of T invariably reduce to one of the special 
cases of (1), apart from the fact that most of them are equivalent [1] up to a given 
level of accuracy. 



Third, although it is clear that T provides a nontrivial contribution to the effec- 



lem 




(i) 
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On the other hand, analyses have been performed comparing experimental data 
with theoretical insights provided by various considerations such as, for example, the 
study of the envelope function for electrons in uniform or slowly graded crystals [1]. 
These, however, do not always reveal [25] precise and unambiguous estimates to the 
question of uniqueness of the parameters a, ft, 7. Some of the appropriate single- 
band PDEM Hamiltonians are the ones of BenDaniel and Duke (BDD) [26] (a = 0, 
(3 = -1), Bastard [27] (a = -1, (3 = 0), Zhu and Kroemer (ZK) [28] (a = -\, 
(3 = 0) and the redistributed model [29] (a = 0, /3 = -\). 

With the correspondence p — > — ih-^ in (1), the time-independent Schrodinger 
equation reads 



rn a (x)—m^(x)—m y (x) + m y (x)—m l3 (x)—m a (x) 
ax ax ax ax 

Eif){x) 



+ V(x)\ip(x) 



(2) 



for some given potential V(x). On setting 



m(x) = moM(x), 



(3) 



where M(x) is the dimensionless form of the mass function, along with h = 2mo = 1, 
we can get rid of the ambiguity parameters by transferring them to the effective 
potential energy of the variable-mass system. Thus using the result 



ax ax ax ax 

did M" M' 2 



(4) 



where a prime denotes derivative with respect to x, equation (2) acquires the form 

did 



Hij(x) = 



+ V eS (x) 



iP(x) = EiP(x), 



dx M(x) dx 

in which the effective potential V e s(x) is seen to depend on some mass terms: 

1 M" M' 2 

V cS (x) = V(x) + -(? + 1) — - [a(a + p + + — 



(5) 



(6) 
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In the following, we shall be interested in bound-state eigenvalues E n , n = 0, 1, 
2, . . . , and corresponding wavefunctions ip n (x), n = 0, 1, 2, . . . . 
Let us consider the intertwining relationship 

V H = H lV , (7) 

where Hi has the same kinetic energy term as H and an associated potential 
Vi,es( x )- If the ground-state wavefunction ip of H is annihilated by the op- 
erator f], i.e., rji/jo = 0, the eigenvalues of Hi are E^ n = E n+ i, n = 0, 1, 
2, . . . , with corresponding wavefunctions r]ip n+ i, since from (7) it follows that 
#i(wWi) = V H ^n+i = K+iivtpn+i) for n = 0, 1, 2, ... . 

Choosing a first-derivative intertwining operator rj = A(x)-^ + B{x) in (7), we 
are led to the restrictions 

A(x) = M" 1 / 2 , (8) 
V eS (x) = \ + B 2 -(AB)', (9) 
V 1:eS (x) = V eS + 2AB'-AA", (10) 

where A denotes some integration constant. 

It is instructive to consider the free-particle case V(x) = Vq of (6). On comparing 
with (9), where we choose X = V and 

we get a constraint relation (3 = —2a — 1 and equation (6) is found to be consistent 
with the PDEM Hamiltonians of BDD (a = 0, (3 = -1) and ZK (a = -\, (3 = 0). 
Specifically, V e s and Vi je ff read 

M" M' 2 
V cS (x) = V -a— + a(a + 2) — , (12) 

IWx) = K 0+ ( a + 5 )_ + ( a+ -)(«--)_, (13) 

which are interchanged under the transformation a — > — (a + |). In particular, for 
a = 0, 

K ff (x) = V , y liCff (,) = y + -_--_ F , (14) 



while for a — — \, 

1 M" 3 M' 2 

V ^= V °+2lp-lM^ ViM^) = V , (15) 

suggesting a duality between the BDD and ZK schemes. It is to be stressed that 
the results (8) - (15) are independent of any choice of M(x). 

To proceed further with the effective-mass Schrodinger equation (5) we need to 
have a precise example for M(x). We can take for instance a deformed hyperbolic 
function 

M(x) = sech 2 qx, q > 0, (16) 
which depicts a solitonic profile. Another acceptable form for M(x) is 

M(x)=(l + r -^) 2 , g>0, (17) 

which has yielded interesting connections [19] with the su(l,l) algebra. In (16) or 
(17), q may be treated as a deformation parameter so that when q — > 0, M(x) — > 1 
in both the cases. 

Let us take for concreteness the result (15), which is in conformity with the ZK 
scheme. We can write for the counterpart of (5) for Hi 

where E' ln = E^ — Vq = E n+ i — Vo = E' n+1 and tp n (x) = r)i() n+ i(x) are the eigenfunc- 
tions of Hi. Using the specific example (16) for M(x) and taking (8) and (11) into 
account, we can write these eigenfunctions as <p n (x) = ^(coshgx^n+i)- Inserting 
this expression in (18), integrating once the resulting equation, changing n + 1 into 
n and rewriting the result in terms of Xnif) = cosh qx ip n (x), where t = tanhqx, we 
find 

a-^-^+J^.o. (19) 

Such an equation coincides with the equation for Legendre polynomials Xn(t) = 

P n (t), n = 0, 1, 2, . . . , provided E' n = q 2 n(n + l). 1 From this, it follows that 

1 Actually the general solutions 2-Fi(— f, f + 1; 1; (1 — tanhgx)/2) of Legendre's differential equa- 
tion (19), corresponding to E' v = q 2 v(i> + 1) (y ^ n), being not convergent for x — » — oo, the 
polynomial solutions are the only acceptable ones. 
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apart from some normalization factors, i]) n {x) oc sechgxP n (tanhgx) and <p n (x) oc 
^P„ + i(tanhgx). We therefore get for the first few normalizable solutions of H and 
Hi, 

ipo{ x ) ^ sechgx, E' = 0, 

tpi(x) oc sech qxtanhqx, fo{x) oc sech 2 qx, E[ = E' 10 = 2q 2 , 
■02 (x) oc sech qx(l — § sech 2 qx), fi(x) oc sech 2 gxtanhgx, 
E' 2 = E[ A = 6q 2 , 

ips(x) oc sech qxta,iohqx(l — | sech 2 qx), <f2(x) oc sech 2 qx(l — \ sech 2 qx), 
E' 3 = E' h2 = Ylq\ 

■04 (x) oc sech qx(l — 5 sech 2 qx + y sech 4 gx) , ^(x) oc sech 2 qx tanh gx 

x (1-| sech 2 qx), E' 4 = E' 13 = 20q 2 . (20) 

We conclude that a free particle placed in an appropriate mass background generates 
bound states in a manner as given above. 

3 Coordinate Transformation 

Equation (5) can be interpreted as a coordinate-transformed Schrodinger equation 
with the coordinate dependence arising from a suitable choice of mass function 
M(x). To this end, let us consider a time-independent Schrodinger equation with 
constant mass under the action of a potential U, 

d 2 



+ U(y;a) 



<Ky) = e<f>(y), (21) 



dy 2 

where a denotes collectively a set of coupling parameters, which may be present in 
U, and e is the energy. 
The change of variable 

y = Xz(x) + v, A,i/6R, (22) 
z(x) = j X ^jM(x') dx', (23) 
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transforms equation (21) into 
1 d 1 d 



+ \ 2 U(\z(x) + v] aj 



X(x) = X 2 ex(x), 



(24) 



VM dx VM dx 

where x(x) = (f>{y{x)). 

To compare (24) with (5), one needs to write the term containing the second- 
order derivative in the former to transform in the same way as in the latter. For 
such a purpose, we exploit the property ^Jm^ = T^T/tr + \it^h and substitute 



dx y/M ~ 2 M 3 / 2 



X(x) = M- l / A (x)i)(x) to reset (24) as 

7M' 2 



d_J_d_ U" 



+ \ 2 U(\z(x) + v\ aj 



if)(x) = X 2 e / ip(x). 



(25) 



dx M dx AM 2 16M 3 
Equations (5) and (25) are seen to coincide (with E = A 2 e) provided V cS is 
identified as 



V cS (x) = \ 2 u(\z{x) + u; aj + 



M" 1M 
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AM 2 16M 3 ' 

From (6) and (26), we thus arrive at an important result 

V(x) = Vi (x; a,P) + V 2 (x; a) , 



(26) 



(27) 



where Vi(x;a,/3) and V 2 (x;a) are 



Fi(x;a,/9) 



a(a + (3 + 1) + (3 + 



V 2 (x;a) = \ 2 u[\z(x) + v;a). 
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M' 2 l. na . M" 
M3--4 (2/3+1) M^' 



(28) 
(29) 



We see at once that for any solution (e, 4>(y)) of the constant-mass Schrodinger 
equation, we can derive a corresponding one (E,i(j(x)) of the effective-mass 
Schrodinger equation with E = X 2 e and the potential V(x) given by the set (27) - 
(29). 

4 Analysis of V(x) 

Of the two parts of V(x), one involves the ambiguity parameters of the T operator, 
which we call Vi(x; a, (3), while the other (i.e. V^xja)) is a coordinate-transformed 
piece obtained as a result of coordinate transforming U according to (22) and (23). 
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4.1 The potential Vi(x; ex., f3) 

Unrestricted a, (3 parameters are often of much use to study the general class of 
solutions for the PDEM Schrodinger equation and to understand their role when the 
system undergoes a transition from a smooth potential and mass step to an abrupt 
situation [11]. In spite of this, an interesting question is whether the presence of V\ 
can be eliminated from V for special choices of a, (3 and/or the mass function M(x). 
The answer turns out to be in the affirmative as our analysis below will presently 
confirm. 

Writing V\(x; a, (3) as 

M' 2 M" 

V^x; a, (3) = f(a, (3)— - g((3)—, (30) 

where f(a,(3) = (a + j) 2 + \{a + l)(2/3 + 1), g{(3) = \(20 + 1), we notice that 
V\ vanishes for (3 = — \ and a = — \, corresponding to any smooth mass function 
M(x). The value (3 — — | is consistent with the redistributed model [29], which is 
a special case of (1), as noted earlier. Although the accompanying zero value of a 
is different from — |, the PDEM kinetic energy operators are known [30, 31] to be 
equivalent over a wide range of a values for suitable width of the transition region. 

For (3 7^ — |, the vanishing of V±(x; a, (3) also takes place for the following mass 
dependence 

f^g: M ~ £ = (l - L ) , (31) 

f = g: M ~ e ±kx , k e R + , (32) 

i.e., either for power and inverse power laws or for exponentially rising and falling 
masses. For the kinetic energy operators of BDD [26], Bastard [27], and ZK [28], 
£ turns out to be — |, — |, —4, respectively. In such cases, M(x) is singular for 
x — > 0. In general (i.e., notwithstanding the vanishing of Vi), the case f — g 
is disfavoured by the above schemes because the underlying constraint equation 
( a + l) + \ ( a + |) (2/? + 1) = is not satisfied by the corresponding parameter 
values. However, exponential dependence on the position coordinate of the mass 
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has been considered [22] viable in semiconductor quantum well structures and also 
in problems of exact solvability in the supersymmetric context. 

4.2 The potential V 2 (x;a) 

The comprehensive nature of our scheme enables us to obtain rather straightfor- 
wardly the relevant forms of V 2 should U for the constant-mass Schrodinger equation 
be available. To show how our scheme works in practice, we consider the following 
illustration of the Scarf I potential, which is well known to be ES on the interval 

U(x;A,B) = [A{A - 1) + B 2 ] sec 2 x - B{2 A - 1) sec xt&nx- A 2 , 

0<B<A-1, (33) 
e n = (A + n) 2 -A 2 , n = 0,1,2,.... (34) 

Using (23), (29) and say, for instance, (17), we easily find that in the presence of 
mass deformation the argument x in (33) is modified according to x — > Xz(x) + v, 
z(x) = x + gtan" 1 x. As a result, V 2 (x; A, B) can be written as 

V 2 (x; A, B) = [A(A - 1) + B 2 ] sec 2 [z(x)} - B(2A - 1) sec[z(x)] tan[z(x)} - A 2 , (35) 

where along with v = 0, we have set A = 1 so that the energy levels are the same as 
given by (34), i.e., E n = e n . 

A graphical representation of (35) is shown in figure 1, where we have considered 
the effects of g-deformed mass function (17) for various q values and compared them 
with the undeformed (i.e. q — 0) case (33). It is to be noticed that V 2 (x] A, B) is 
not singular provided 

-ig + ,)<tan-x<ig-x). (36) 

Hence it is defined on the interval (—x q , x q ), where x q G (o, |) is the solution of the 
transcendental equation tan -1 :^ = ^ (f — As q increases, x q decreases from 
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xq — | to = 0. One finds, for instance, xi/io — 1.47335, X1/2 — 1.14446, and 
xi ~ 0.860334 for the curves displayed in figure 1. 

A similar procedure can be used by taking as inputs the ES potentials known in 
quantum mechanics [32, 33] or generalizations thereof. In this way, we can get in a 
unified way many of the exact results that have been derived elsewhere by some ad 
hoc techniques [10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22]. 

As two new examples, let us mention 

• Kratzer potential 



V 2 (x; 1 ) = 4 7 2 --r^-^T, 7>0, 
16z z (x) z(x) 

E n = 4 7 2 - I— , n = l,3,5,.. 



(2n+l) 2 ' 



Symmetric QES potential 



V 2 (x;A) = -A 2 smh 2 [z(x)]- Acosh[z(x)] + -A + -, A > 0, 
4 2 4 

E = 0, E 1 = A, 



where we have taken v = for simplicity and we have fixed A = 1 without any loss 
of generality. 

We remark that one of the couplings in Kratzer is fixed at — ^ and so not free. 
This entitles it to be classified as a CES potential. Further we observe that the run 
of n is restricted to odd values only because physically acceptable wavefunctions 
then have the correct x 3 / 4 behaviour at the origin (for a discussion of wavefunction 
behaviour for strongly singular potentials see [34, 35] and references quoted therein). 
The symmetric QES potential [36] , on the other hand, is a special case of Razavy 
potential [37] with two known eigenstates. It is needless to say that broader classes 
of potentials, such as the ones obtained in [38, 39, 40, 41, 42, 43, 44], may also be 
mass deformed following the approach prescribed in this paper. 
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5 Conclusion 

We studied in this Letter a PDEM quantum Hamiltonian in one dimension guided by 
the kinetic energy operator T of von Roos. The free-particle problem is analyzed first 
in the presence of a sech 2 -mass background and generation of bound states is noted 
by exploiting a first-order intertwining relationship. The non-trivial nature of this 
result should be stressed. In a second step, the accompanying potential V(x) of the 
PDEM Hamiltonian is found to reduce to two terms — one involving the ambiguity 
parameters of T and the other emerging from the coordinate transformation of 
the constant-mass Schrodinger equation. The advantage of our scheme is that, for 
a given mass function, a knowledge of some constant-mass Schrodinger potential 
allows a full access to the associated potential in the effective-mass Hamiltonian. 
This provides a unified treatment of all mass-deformed potentials corresponding to 
ES, QES or CES potentials known in the constant-mass case. We illustrated our 
results by an appropriate choice of the mass function as applied to the well-known 
Scarf I potential. 
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Figure Caption 

Fig. 1. Comparison between the potential V 2 (x; 3, 1) of Eq. (35) for q = 0.1 (dashed 
line), q = 0.5 (dotted line), or q = 1 (dot-dashed line) and the Scarf I potential 
U(x; 3, 1) (solid line). 
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